WHAT IS CLAIMED IS : 

1. A method for retrieving a data image similar to a query image 
using a feature distance calculated by combining one or more color distances 
and one or more texture distances by considering human visual perception 
attributes. 

2. A method for retrieving a data image similar to a query image 
in an image database containing a plurality of data images, the method 
comprising: 

(a) calculating a plurality of color distances and a plurality of texture 
distances between a query image and each data image in the image database; 

(b) weighting each of the calculated color distances and texture 
distances with a respective predetermined weighting factor; 

(c) calculating a feature distance between the query image and each 
data image by combining the weighted color distances and the weighted 
texture distances by considering human visual perception attributes; and 

(d) determining the data image similar to the query image using the 
feature distance. 

3. The method of claim 2, before step (a), further comprising: 

(pa-1) segmenting the query image and each data image into a plurality 
of first regions using a plurality of color features; and 

(pa-2) determining a plurality of sample regions in. the query image 
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and each data image for extraction of a plurality of texture features. 

4. The method of claim 3, wherein step (a) comprises: 

(al) generating a plurality of color vectors of the first regions using the 
color features and calculating a plurality of color distances; and 

(a2) generating a plurality of texture vectors of the sample regions 
using the texture features and calculating a plurality of texture distances. 

5. The method of claim 3, wherein step (pa-1) comprises 
quantizing a plurality of color vectors of the query image and each data image. 

6. The method of claim 5, wherein quantizing the color vectors 
comprises: 

(pa- 1-1) performing a peer group filtering on the query image and each 
data image for noise removal and smoothing effects; and 

(pa- 1-2) clustering a plurality of filtered pixel values of the query 
image and each data image using a generalized Lloyd algorithm. 

7. The method of claim 5, further comprising: 

defining a J-value indicating a color uniformity in each pixel of a 
plurality of pixels of the query image and each data image, which have 
undergone quantization; 

storing the J-value in each pixel of the query image and each data 
image to obtain a plurality of J-images; 
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segmenting each J-image into a plurality of second regions by a 
predetermined segmentation method; 

repeating the segmentation of each J-image to obtain a map of one or 
10 more over-segmented regions for each J-image; and 

obtaining a final map for each J-image by merging the over-segmented 
regions based on a correlation of color. 

8. The method of claim 7, further comprising indexing a feature 
vector space by a representative color and a percentage of the representative 
color in each second region. 

9. The method of claim 6, wherein step (pa- 1-2) comprises 
applying a predetermined algorithm to increase the number of resulting 
clusters or to merge the resulting clusters. 

10. The method of claim 4, wherein the color features are 
expressed by a color feature descriptor f c (I k ) with a representative color 
vector and a percentage of the representative color vector for each first region. 

11. The method of claim 10, wherein the color feature descriptor 
f c (I k ) is expressed by: 

fc(h> {(^1 > Ph )> (^2 > Pk2% - > (fiw k > PkN k )} 

wherein k is a positive integer indicating a serial number of each region, 
5 is an i-th representative color vector of a k-th region (i = 1,2,.., N), P h is a 
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percentage of the i-th color representative color vector in the k-th region, and 
N k is the number of the representative color vectors in the k-th region. 

12. The method of claim 8, wherein indexing the feature vector 
space comprises: 

assigning a plurality of representative colors to a plurality of grid 
points in a color space having a grid structure; and 

storing the result of the assignment as a table in a database. 

13. The method of claim 11, wherein the color distance between 
each color vector of each region of the query image and each color vector of 
each region of each data image is calculated using: 

^ n 2 ^ N 2 

d c uM = £ pit 2a ^j PuPzj 

1=1 1=1 1=1 J=l 

wherein when T d is a maximum distance by which a similarity of two colors 
is determined, a is a predetermined scaling coefficient, d max = a T d , and d y is 

the Euclidean distance ||c, - c y || between two color vectors c J and c 3 , such 
d 

that a v = 1- -7 1 - if d y < T d , and a y = 0 if d y > T d . 

max 

14. The method of claim 4, wherein step (a2) uses a Gabor 
function. 

15. The method of claim 14, wherein, in step (a2), the texture 
vectors of the plurality of sample regions are generated using the Gabor 
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function having N frequency channels and M orientation channels, where N 
and M are both predetermined positive integers. 

16. The method of claim 15, wherein the texture features are 
expressed by a texture feature descriptor f t (I k ) : 

/,(/*)= {(rn kn9 a kn ),(™ kn ,<T kl2 ), ...,(iw«A/^HA/X( w *2i^it2i)> -'("V —>( m *NM> a *NM)} 

wherein m hj is a mean value of a plurality of pixel values of the i-th frequency 
channel and the j-th orientation channel for a sample region I k , and a hj is a 
deviation of the pixel values of the i-th frequency channel and the j-th 
orientation channel for the sample region I k . 

17. The method of claim 16, wherein, in step (a2), the texture 
distance between each color vector of each region of the query image and each 
color vector of each sample region of each data image is calculated using: 



N M 



;=1 7=1 



N M 

+ IS 

1=1 j=l 



and 



18. The method of claim 3, wherein step (pa-2) comprises: 

(pb-1) obtaining a grid map of the query image and each data image; 



(pb-2) obtaining a texture sample of a desired size for each sample 
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region based on the grid map. 

19. The method of claim 18, wherein step (pb-1) comprises: 

(pb-1-1) taking a rectangle M(iJ) having largest side lengths for a 
sample region, wherein 0 <i <n and 0 <j <m; 

(pb-1 -2) dividing the rectangle M(i 9 j) into a plurality of sub- 
rectangles each having a / x / size; and 

(pb-1 -3) obtaining a grid map for the sample region using a 
predetermined function which outputs 1 if all the divided grids belong to the 
sample region, and otherwise, outputs 0, the predetermined function expressed 
by: 



\\ if p y (x,y) c I k for\f(x,y) 



0 otherwise 



wherein 0 < i < 



1, 0<7< 



1, 0<x</-l, 0<y<l-l, and 



p y (x,y) is a point that belongs to rectangle M(il + x, jl + y) . 

20. The method of claim 19, after step (pb-1 -3), further 
comprising: 

(pb-1 -4) transforming the grid map obtained in step (pb-1 -3) to a grid 
distance map by repeating computations with: 
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{M^iUj)}^ = {M% ld {Uj)} n ,M d grid {iJ) = {M^(/,y)}„, an d 

(pb-1-5) performing region growing in every direction from a seed of 
10 the grid map, which has the largest distance in the grid and is expressed by 
(a,Z>)= « {M^(/j)}, to satisfy the relation M d grid (U)>0 until a 
maximum area rectangle is extracted. 

^ 21. The method of claim 20, further comprising fitting the 

g maximum area rectangle extracted in step (pb-1-5) to the desired size of the 

.ass, 

jfU texture sample by wrapping. 

JL, 22. The method of claim 20, further comprising fitting the 

q maximum area rectangle extracted in step (pb-1-5) to the desired size of the 

■yy 

O texture sample by mirroring. 

23. The method of claim 3, wherein step (b) comprises: 

(b-1) placing each of the color distances and the texture distances in a 
2-dimensional vector space, each vector space defined by the respective 
distances and associated predetermined weighting factors; and 

5 (b-2) projecting the result of the placement onto the 2-dimensional 

vector spaces onto a 1 -dimensional distance space based on the human visual 
perception mechanism. 
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24. The method of claim 2, wherein the predetermined weighting 
factor to the color distance is determined based on a distribution of 
representative colors. 

25. The method of claim 23 , wherein the predetermined weighting 
factor to the color distance is determined by: 

N i 

®c = l-Z^ lo Sio(~) 

i=l Pqi 

wherein p is a percentage of an i-th representative color of one of the first 
regions of the query image. 

26. The method of claim 23 , wherein each texture distance is defined by: 

dl{I q J x )=a A ^ d t {I q J x ) 

wherein I q denotes the query image or each first region of the query image, s 
denotes a sample region of a desired size, AQ) denotes the area of the sample 
region of the desired size, countQ) is the number of wrappings done to obtain 
the sample region of the desired size, and a is a constant. 

27. The method of claim 23, wherein the predetermined weighting 
factor applied to each texture distance is determined based on an area of an 
initial sample region extracted from the query image and the area of a sample 
region extracted from each data image. 
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28. The method of claim 27, wherein the predetermined weighting 
factor applied to each texture distance is determined by: 

1 

a q 

wherein I q denotes the query image or each first region of the query image, s 
denotes a sample region of a desired size, A(>) denotes the area of the sample 
region of the desired size, count (•) is the number of wrappings done to obtain 
the sample region of the desired size, and a is a constant. 

29. The method of claim 23, before step (b-2), further comprising 
normalizing each of the color distances and the texture distances. 

30. The method of claim 29, wherein normalizing each of the color 
distances and the texture distances is performed using a Gaussian 
normalization. 

31. The method of claim 30, wherein normalizing each of the color 
distances and the texture distances comprises: 

wherein v k is a Gaussian sequence, performing a normalization within 

v rn k - jU k 

a range of [-1, 1] using v mk = — ~ based on a mean value //^and a 

3<r k 

deviation a k of the sequence v k ; and 

mapping the result of the normalization into a range of [0, 1] using 
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32. The method of claim 31, wherein normalizing each of the color 
distances and the texture distances comprises normalizing the texture distances 
by updating the mean value ju t and the deviation a t of the texture distances 

by excluding a largest texture distance until the condition of kx— < — , 

wherein £ is a constant, is satisfied. 

33. The method of claim 23, wherein step (b-2) comprises 
projecting the color distances and the texture distances onto a 1 -dimensional 
distance space using: 



d(I q J l )=W c d c (IJ } ) 



( 2 _,d< t (I q jS" 
1+ —tan 



vs in t 
wherein W c = — , and W t = . 



+ W t d< t <J q ,I,) 



( 2 _, d c {I q J^ 
1 + — tan 1 



V * d' t {I ,1^) 
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